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Abstract. - We investigate the structure of colloid-polymer mixtures by calculating the struc- 
ture factors for the Asakura-Oosawa model in the PY approximation. We discuss the role of 
potential range, polymer concentration and polymer-polymer interactions on the colloid-colloid 
structure. Our results compare reasonably well with the recent experiments of Moussaid et. 
al.^ for small wavenumber k, but we find that the Hansen- Verlet freezing criterion is violated 
when the liquid phase becomes marginal. 



Suspensions of sterically stabilized colloidal particles in an organic solvent containing non- 
adsorbing polymer have been widely studied experimentally and theoretically because of their 
remarkable phase behaviour. The non-adsorbing polymer coils induce an attractive depletion 
interaction between the colloidal particles, the range and depth of which can be easily "tuned" 
by varying the polymer molecular weight (and hence its radius of gyration Rg) and the polymer 
concentration rip. This makes polymer-colloid mixtures an important experimental model 
system to study the role of the range of attractive interactions on phase behaviour. 

If Rc is the radius of a spherical colloidal particle, and ^ = Rg/Rc denotes the size ratio, it 
was predicted theoretically |l|, |^, and shown experimentally (|] and by simulations Q] that the 
phase diagram in the (rjc, rip) plane (where rjc = AirricRc/S denotes the colloid packing fraction) 
exhibits colloid "gas", "liquid" and "solid" phases, when ^ >~ 0.3, while the intermediate 
"liquid" phase becomes "marginal" for f 0.3, and disappears completely for smaller size 
ratio^(^. The previous theoretical work mostly focuses on calculating free energies and the 
resulting phase diagrams and no attempts were made to investigate the structure of the colloid- 
polymer mixtures, as embodied in the partial structure factors Sa/3{k) (1 < Q < /3 < 2), where 
the indices 1 and 2 refer to the colloid and polymer species respectively). This is because the 
van-der-Waals like mean-field theories in[^ |^ do not readily yield structural information, and 
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(1) There is stiU some discrepancy between theory and simulations which predict a marginal liquid 
for 5 ~ 0.4|^, ^ , and experiment, which finds a marginal liquid at ^ ~ 0.24j^. 
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because until very recently no experimental data were available. This gap has been filled by the 
two-colour dynamic light-scattering measurements of Moussaid et. al. ^ , which provide data 
for the colloid-colloid structure factor, S{k) — Sii{k), over a limited range of wavenumbers fc, 
for thermodynamic states at triple-point conditions, and for three values of the size ratio ^. In 
particular, they concluded that under triple point conditions the amplitude of the first peak 
of the structure factor was always near 2.5 while the small k behaviour showed a steep rise 
reminiscent of critical behaviour. In other words the local fluid structure at the triple point 
is insensitive to the potential range, while the large scale structure is sensitive to this range. 
These finding are analyzed here within a simple colloid-polymer model. 

In this letter we apply the familiar Percus-Yevick (PY) closureQ to the Asakura-Oosawa 
(AO)|0] model of a colloid-polymer mixture to extract the partial pair distribution functions 
gapif), and corresponding structure factors. The AO model treats a colloid-polymer system 



as a non- additive mixture of spheres of radius Rn = Rc (the colloids) and point particles 



(the weakly interacting polymer coils) which are excluded from a sphere of radius R12 — 



Rc + Rg around each colloid Q_ . This binary system is characterized by three total and direct 
correlation functions, haffir) = gapir) — 1, and Ca/sir), the Fourier transforms of which satisfy 
three coupled Ornstein-Zernike relationsra: 



{k)h^p{k). 



We have solved eqs (^ numerically, subject to the PY closure relations: 

Cap{r)^ 0,r>Rai3. (2) 

Note that (||) implies C22{r) — 0, since R22 — (point particles). The 2nd and 3rd OZ relations 
(0) accordingly simplify to yield the relations: 

S'i2(fc) = ^Jnrmhi2{k) = y/nin2Ci2{k)Sii{k) (3) 
S22{k)= l + n2h22{k) =l + nin2cl^{k)Sn{k). (4) 

Note that 5*22 (fc) 7^ 1 (and hence h22{r) ^ 0), despite the fact that the polymer coils do not 
interact (point particles) , due to correlations induced by the presence of the colloidal particles. 
Thus (722 (^) provides a measure of the "void structure" of the annealed, porous "medium" 
created by the colloidal spheres. 

Contrary to the case of additive hard-sphere mixtures[^, the coupled OZ/PY equations 
admit no analytic solution in the present non- additive caseQ; an analytic solution for the 
AO model only exists in the physically irrelevant case of negative non-additivity (which would 
amount to Rg < 0!) Accurate numerical solutions are, however, readily obtained and 

results are reported below, and compared to the experimental data and to the PY results 
for the effective one-component system of colloidal particles interacting via the AO depletion 
potential^]. The effective depletion interaction follows from the gain in free volume, accessible 
to the non-interacting polymer coils, due to overlaps of exclusion spheres of radius R12 for 
any configuration of the non-overlapping colloid spheres of radius Rn < Ri2- The AO pair 
potential follows from a second-virial-like approximation whereby only pair-wise overlaps are 

(^) By non-additive we mean that the R12 ^ |(^ii + -R22). 

{^) Since interactions between polymer- coils are neglected, their activity ap — exp{l3fj,p} is propor- 
tional to the effective polymer concentration — Np/V' , where V' is the free volume accessible to 
the polymer coil; V' is a fraction aV of the total volume due to the presence of the colloidal particles 
from which the polymer coils are excluded. 
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considered, thus neglecting more than two-body effective interactions resulting from overlap 
of three or more exclusion spheres. This is strictly valid only for ^ < S,o — 2/\/3 — 1. For 
^ > ^0 the two-component AO model can no longer be exactly mapped onto an effective 
one-component system of colloidal particles with strictly pairwise interactions. 

For a one-component system the long wavelength (fc 0) of the structure factor S{k) is 
given by the familiar relation[^: 

Vim^S{k) = ricksTxT- (5) 

where xt is the isothermal compressibility of the effective one-component system. The 
situation is more complicated in the two-component case, where the (fc 0) limit of the 
Capik) and of the Sap{k) are governed by the derivatives of the chemical potentials of the two 
species, fiap — nfj{d{fia/ ksT) / dni3)T,na ■ For the AO model in the PY approximation, where 
£22(0) = 1 — /i22 = 0, implying ^22 = 1, one finds in particular that: 

\imSn{k) = ^ , (6) 

k^o fill - M12M21 

whereas (with n ^ ni + n2 and = na/n) 

1 

nksTxT = lini 



fe-^o 1 — Xincii{k) — 2xiX2'nci2{k) — x'^nc22(k) 

a;2 + 4- 2^12) ' '^'^^ 

showing that the long wavelength limit of the colloid-colloid structure factor is not directly 
proportional to the osmotic compressibility of the colloid-polymer mixture. Note that the PY 
closure implies that /X22 = 1, which is tantamount to assuming that the free volume V is 
unaffected by the presence of the non- interacting polymer-coils, i.e. is independent of n2] this 
is precisely the basic assumption underlying the mean field theories of the phase diagram||^, H- 
The overall behaviour of the three partial structure factors is illustrated in fig. |], where 
a comparison is made with the analytic PY results for an additive mixture of spheres and 
point-particles (corresponding to the Rg = limit of the AO model; Sii{k) is then unaffected 
by the point particles and thus equivalent to a one-component hard-sphere model). The 
non-additivity is seen to strongly affect Si2{k) and S22{k), while the colloid-colloid structure 
factor, S'ii(A;), is less sensitive. The greater sensitivity of 5*12 (fc) and S22{k) may be traced 
back to the significant change in void structure when the excluded volume radius R12 increases 
(non-additivity) . 

Fig. ^ compares S'ii(fc) calculated for the AO model to the experimental data of Moussaid 
et.aZ.^ for the three size-ratio's they used: ^ — 0.57, which corresponds to a situation where 
the critical and triple points are well separated, ^ = 0.37, where the critical point and triple 
point almost merge and the liquid is "marginal" , and ^ = 0.24, where there is no stable 
liquid region and the critical point becomes metastable w.r.t. the melting line (Note that 
in the experiment, ^ = 0.24 (and not ^ = 0.37) corresponds to the marginal liquid.). The 
colloidal packing fractions, rjc, and size-ratios, ^, were taken from the published experimental 
values while the polymer packing fractions, rjp = AnnpRg/Z, (which depend on Up and ^) were 
adjusted to fit the experimental structure factors, Sii{k), at k = 0.7/RjPJ. Keeping in mind 



('') We note that taking the experimental polymer number density as fixed and varying the effective 
size-ratio parameter, ^, would be more consistent with fitting to experiment since it is not clear that 
^ should be exactly equal to Rg/Rc. For example we find we can fit the small k behaviour of S'ii(fc) 
at the experimental polymer density for the shortest polymers by changing ^ from 0.24 to 0.28. This 
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Fig. 1. - Comparison of the structure factors of the 2-component AO model and an additive mixture 
of hard spheres and point particles, calculated in the PY approximation for rjc — 0.333, rjp — 0.1 and 
^ = 0.24. 

Fig. 2. - Comparison of the experimental colloid-colloid structure factors of Mousaid et. al. jsj with the 
AO model in the PY approximation. The values of the colloid packing fraction, rja, and the size ratio 
^ are taken from experiment. The polymer packing fraction, rjp, is fit to the experimental value of 
SwikRc) at the smallest accessible wavenumber (kRc = 0.7). The resulting polymer packing fractions 
are rjp = 0.225,0.15 and 0.089 respectively, compared to the measured values of r]p — 0.13,0.13 and 
0.1 respectively. The 2nd theoretical line at ^ = 0.24 demonstrates the insensitivity of the 1st peak 
versus the sensitivity of the low k behaviour of the colloid-colloid structure factor, 5*11 (fc), to changes 
in r]p. 



the fact that PY theory tends to overestimate the amplitude of the main peak, the agreement 
for ^ = 0.57 is satisfactory, but for ^ = 0.37 and ^ = 0.24 the amplitude of the main peak 
drops considerably, in contrast to the experiment where the main peak is insensitive to the 
size ratio ^. 

For a given colloidal packing fraction, rjc, the amplitude of the first peak of the colloid-colloid 
structure factor is rather insensitive to varying rjp (or ^). While shortening the range of 
the attractive interactions widens the melting transition, bringing the liquidus line to lower 
colloidal packing fractions, we find that the amplitude of the first peak of the colloid-colloid 
structure factor is determined primarily by the colloid packing fraction, ijc, and not by the 
proximity of the freezing transition. In contrast, the small k behaviour of the structure 
factor is strongly affected by the proximity of a critical point or a spinodal line. These 
conclusions are consistent with recent simulations of one-component systems with short-range 
potentials fl^, [l^ , but inconsistent with the experiments §. 

Perhaps the most notable approximation in the AO model is the complete neglect of 
polymer-polymer interactions To obtain a rough estimate of the influence of such 

interactions on the pair structure of a colloid-polymer mixture we use the following model for 
an effective polymer-polymer pair potential: 

V22{r)^eeM-r/2Rg)^, (8) 



seems to be consistent with recent findings on the interactions between a gaussian polymer and a hard 
wallO, but is inconsistent with another computer simulation study of the AO modelM. 



Structure of Colloid-Polymer Mixtures 



5 




Fig. 3. - Pair distribution functions for a colloid-polymer mixture with repulsive polymer-polymer 
interactions of the form of equation (^) with r]c = 0.333, rjp — 0.13 and £, — 0.24. While the 
polymer-polymer radial distribution function (722(?') is strongly affected, the colloid-colloid radial 
distribution function (?ii(r) is very weakly affected by interactions. 

Fig. 4. - Comparison of the colloid-colloid structure factor for the AO model and the AO potential 
for ^ = 0.37 , rjc — 0.404 and different values of the polymer packing fraction rjp. The AO potential 
polymer fugacities are converted to polymer packing fractions through the scaled-particle free- volume 
factor a{ric, 



where r is the distance between the centers of mass of the two polymer coils, and e is an energy 
scale which is expected to be of order fesT, since the effective interaction is largely of entropic 
origin. This form of the effective potential is consistent with the findings of direct simulations 
of hard-sphere polymer chains [p^. Eq (|^) is a simple parameterization of the free energy cost 
linked to the interpenetration of two polymer coils; the range of such an effective interaction 
is obviously of the order of the radius of gyration, and depending on the solvent-monomer and 
monomer-monomer interactions, e could be positive or negative. In a good solvent e is expected 
to be positive (effective repulsion). It is easy to incorporate the pair potential (||) in the PY 
theory eq. (H). The effect of turning on the polymer-polymer interactions is illustrated by the 
radial-distribution functions in fig. ^ The interaction (^, turns out to have a weak infiuence 
on the colloid-colloid structure, while, as expected, the polymer-polymer correlations, which 
are only indirect in the point particle case, are strongly affected by V22ir). Thus it appears 
that the discrepancy between the AO model and the experiments does not lie in the neglect 
of polymer-polymer interactions. 

A final question to address is the role of many-body interactions, naturally taken into 
account in the two-component system, but neglected in the popular one-component approx- 
imation to the AO model: the hard-sphere system with an attractive AO potential. The 
comparison is somewhat complicated by thermodynamics, since the strength of the AO 
potential is directly proportional to the polymer activity Op, while the AO model is most 
naturally interpreted in terms of the polymer density Up (or polymer packing fraction r/p). 
Since the polymers are non-interacting, however, the two can be directly related through the 
free-volume fraction V' /V = a. In fig. IS we compare the PY results for the 2-component 
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AO model in the marginal liquid state — 0.37) to the 1-component AO potential with 
the rjp scaled by the appropriate free-volume fraction a for which wc take the scaled-particle 
approximation, which has been shown to be sufficiently accurate for the qualitative behaviour 



we are addressing Q. Significant differences between the structure factors calculated within 
the one and two-component models are observed for identical values of rjp, particularly at small 
k. Use of the AO potential will lead to a severe underestimate of the polymer density, Up] the 
discrepancy becomes progressively worse as ^ increases. For example, for ^ = 0.57 the low-k 
fit to the experimental S{k) yields rjp = 0.088 for the AO model but rjp 0.01 for the AO 
potential. We note that the rise in the structure factor of the AO potential system is consistent 
with the location of the spinodal line calculated in the mean- field theories |Q ^, while the rise 
in the structure factor of the AO model is more consistent with the experiments. Since the 
same PY level of approximation is used to calculate both structure factors, the disagreement 
suggests the inadequacy of the AO potential or the assumption of pair-wise additivity in the 
reduction to an effective one-component system. It must be stressed however, that the PY 
closure is far from exact, and may be inadequate to cope with a short-range attractive well, as 
in the AO potential, particularly for the smaller values of ^. On the other hand PY is expected 
to be more accurate for purely repulsive interactions as in the 2-component AO model. We 
are presently examining this point using thermodynamically self-consistent closures jl^. 

As pointed out earlier, the k limits of the one-component colloid-colloid structure 
factor S{k), and its two-component equivalent Sii{k) do not have the same thermodynamic 
meaning. For ijc = 0.404, rjp — 0.15 and ^ = 0.37, the osmotic compressibility derived from 
eq. (||) for the effective one-component system {S{k) = Sii{k)), and from eq. (|^) for the 
original two-component AO model are ncksTxT — 0.49 and ncksTxT = 0.035 respectively. 

We may draw the following main conclusions from the present work: 

a) In contrast to the experiments of Moussaid et. al. ||] we find (both with the 2-componcnt 
AO model and the 1-component AO potential) that the first peak of the colloid-colloid structure 
factor and the related local fluid order is primarily determined by the colloid volume fraction, 
and not by the the proximity to the melting transition. Thus the Hansen- Verlet criterion, 
which states that the first peak of the structure factor at freezing is near 2.8, is violated as 
the range of the potential is shortened to the point where the triple point vanishes. This 
discrepancy with experiment may be related to novel effects arising from an interplay between 
the triple and critical points, which for ^ = 0.24 are in close proximity (marginal liquid). 

b) In agreement with the experiments^, we find that the large scale fluid structure and 
related small k limit of the structure factor is strongly affected by the proximity of a spinodal 
line or a critical point, and thus depends significantly on the polymer density. The small k 
limit of the structure factors is in fact well described by the AO model. 

c) Polymer-polymer interactions appear to have little impact on the colloid-colloid pair 
structure, at least within the naive ansatz of eq. (^. 

d) The one-component description, based on the effective AO pair potential, appears to be 
insufficient to yield a quantitatively accurate description of the colloid-colloid structure factor, 
particularly at small k. It is consistent with the phase-diagrams calculated in the mean-field 
theories, and similarly underestimates the polymer packing fraction for larger size ratios. 

To ensure consistency of the theory, we are presently extending the PY work to other, 
thermodynamically consistent, integral equations which we will use to determine the phase 
diagram of the two-component model. We are also seeking a first principles approach to the 
effective polymer-polymer and polymer-colloid interactions. 

(^) The overall free- volume is well represented by the scaled-particle theory expression for awt, but 
the thermodynamics of theories based on free volume, such as those in ref [^, depend critically on 
d^a/dnl, for which the scaled-particle approximation is less accurate p^. 
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